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Chiral phase transition is investigated in an SU (3 )l x SU (3)i^ symmet¬ 
ric vector meson extended linear sigma model with additional constituent 
quarks and Polyakov loops (extended Polyakov quark meson model). The 
parameterization of the Lagrangian is done at zero temperature in a hybrid 
approach, where the mesons are treated at tree-level, while the constituent 
quarks at 1-loop level. The temperature and baryochemical potential de¬ 
pendence of the two assumed scalar condensates are calculated from the 
hybrid 1-loop level equations of states. The order of the phase transition 
along the T = 0 and = 0 axes are determined for various parameteriza¬ 
tion scenarios. We find that in order to have a first order phase transition 
at T = 0 as a function of a light isoscalar particle is needed. 

PACS numbers: 12.39.Fe, 12.40.Yx, 14.40.Be, 14.40.Df, 14.65.Bt, 25.75.Nq 


1. Introduction 

In [1] it was shown through a zero temperature analysis that qq scalar 
states, such as the oq, and the two /o’s are preferred to have masses 
above 1 GeV. Similar results was obtained with qq states in [2], while in 
[21011] it was shown by using tetraquarks instead of qq states that the 
(tetraquark) scalar masses are in the range 0.6 —1.0 GeV. These results sug¬ 
gests that the physical states ao(1450), iCQ(1430), /o(1370), and /o(1710) 
(or /o(1500)) are predominantly qq states, while ao(980), iCQ(800), /o(500), 
and /o(980) are predominantly tetraquark states. However, states with the 
same quantum numbers do mix, thus the physical scalar particles are mix¬ 
tures of qq and tetraquarks states. In the current case of the extended linear 
(T model (ELcrM) we have only one scalar nonet, thus we can describe one 
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qq, one Kq and two /o (which will be denoted by /q ' ) particles. Conse¬ 
quently, one of the most interesting question is that which physical states 
our fields predominantly are if we investigate the finite temperature/density 
behavior of our model additionally to the zero temperature properties. 

The lightest isoscalar = O"'"'') qq state, /o (also called a) is strongly 
related to the non-strange condensate. Since the larger is the sigma mass 
compared to the mass of its chiral partner (the pion) the larger is the tem¬ 
perature at which m/g approaches m„ in the chiral symmetry restoration, 
one would expect that a large nif^ mass results in a large pseudocritical 
temperature (Tc) at zero baryochemical potentials. On the other hand it 
is a common expectation that the chiral phase transition is of first order 
as a function the baryochemical potential (/r^) at T = 0, and since with 
increasing mjg mass the transition weakens, at some point it is possible 
that the transition becomes crossover This suggest that for a good 

thermodynamic description a small ruf^ mass is needed, and indeed as it 
turns out our approach supports this requirement. 

The paper is organized as follows. In the next subsection we intro¬ 
duce the model by giving the Lagrangian. In Sec. the determination of 
the model parameters is shown, while the description of the approximation 
used to calculate the grand potential together with the field equations are 
presented in Sec. 


2. The Model 

The Lagrangian we shall use has the following form: 

L = Tr[(D^M)+(D^M)] - mgTr(M^M) - \i[Ti{M^M)f - X2Tv{M^Mf 
+ ci(det M + det M^) + Ty[H{M + M^)] - ^ Tr(L2 + rI^) 


+ Tr 

/ii 


^ + ^]{Ll + Rl 


■92, 




+ -j- Tr(MtM) Ti{Ll + rD + h2 Ti[{L^Mf + {MR^f] 
+ 2/i3 Tr(L^M/?^M^) -b ^ - M] T 


( 1 ) 


The covariant derivatives above are given by 


D>^M = df^M-igi{Lf^M-MR^^)-ieA^[T3, M], (2) 

where = ggG^Ti, with Tj = Ai/2 (i = 1,..., 8) denoting the SU{3) group 
generators given in terms of the Gell-Mann matrices A*. Here M = Mg -b 
Mps stands for the scalar - pseudoscalar fields, R^ = V^^—A^ 
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for the left and right handed vector fields (which contain the nonets of vector 
{Va) and axial vector (^a) meson fields), Ae is the electromagnetic held, 
while Gf are the gluon helds. The held strength tensors are {Q G {L, R}) 

QPU ^ QPQU _ - ieA'^ATz, Q^]} , (3) 

while the external helds related to the scalar and vector helds are H = 
^dmg{hoN, hQ]\f, y/2hos), ^ = diag{6N, Sn, ^s) (For more details on the 
model see pQ). 


3. Setting the Lagrange parameters 

There are 15 unknown parameters in Eq. which will be present 
in the held equations at hnite temperature and/or density, namely, mo, 
Ai, A 2 , Cl, mi, hi, /i 2 , hs, Sn, ^s, (i>N, (ps, 9 f, 91, 92 - From this set Sn 
can be melted into mi thus leaving 14 unknowns. These parameters are 
determined similarly as in [T], that is we calculate values of various masses 
and decay widths at tree-level and compare them with the corresponding 
experimental value taken from the PDG [7j through the minimalization 
process of Ref. [8]. It is important to note that we artihcially increased the 
errors of the PDG to a minimum level of 5%, since we do not expect that 
our model to be more precise. Another important points are that now we 
use a different anomaly term in Eq. ([^ (the term proportional to ci) as 
was presented in [1], we ht the total width in case of ao(980) instead of the 
amplitudes, we include the /o masses and decay width into the global ht, we 
take into account of the effects of the fermion vacuum huctuations, case of 
which the expression of the (pseudo)scalar masses are modihed. Moreover, 
since now we also included the constituent quarks in the isospin symmetric 
limit, we use two additional equations to their tree-level masses with the 
values mu,d = 330 MeV, and nis = 500 MeV. 

The scalar meson sector below 2 GeV contains more physical particles 
than we can place into one qq nonet (consisting of oq, Kq, Jq, f^), since 
in nature two oq, two Kq and hve /o particles exist in that energy range. 

These particles are the ao(980) and ao(1450) (denoted by a^^), the Kq{800) 

and it'Q(1430), (denoted by the /o(500) (or cr), /o(980), /o(1370), 

/o(1500) and /o(1710), (denoted by /q'"^). Accordingly there are 40 par¬ 
ticle assignment possibilities to pair the physical particles to the members 
of the nonet in the model. We performed a fit for all the assignments 
and ordered them according to their x^ values. The results of the 5 best 
solution along with the particle assignments in two cases (with and without 
the fermionic vacuum fluctuation) are shown in Table By a similar fit¬ 
ting procedure in [T] we argued that the best assignment without fitting the 
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particle assignment 

x"" 

^red 

particle assignment 

x"" 

^red 


45.0 

3.0 


46.2 

3.1 


51.7 

3.4 


52.7 

3.5 


51.7 

3.4 


54.1 

3.6 

«o^o^/o/o^ 

60.4 

4.0 


60.7 

4.0 


61.3 

4.1 


61.7 

4.1 


Table 1. and Xred = X^/-^dof values (A^dof = 15) for the first five best solutions 
of the fit together with the particle assignment without (left part) and with (right 
part) the fermionic vacuum fluctuations. 


is the while we reasoned that should correspond to 

Now it seems that the situation changes since Table one can see that the 
two best assignment in case with and without using the fermion vacuum 
fluctuation are agi^Q^/g/g and agi^g^/g/g, respectively. This suggest that 
using only zero temperature quantities in this model is not enough to point 
out uniquely a single particle assignment. Thus we investigated the proper¬ 
ties of the different particle assignments at finite temperature/densities as 
well. 


4. Finite temperature field equations 

In our model there are four order parameters, which are the two chiral 
condensates (j)^ (non-strange) and (ps (strange) and the two Polyakov loop 
variables 4? and $ (For the introduction of the Polyakov loop variables and 
their potential see El). The held equations, which determine the dependence 
on T and hb = 3/iq of the order parameters are given by the minimalization 
of the grand canonical potential (see eg. m), 

d(j)N dips 94 > ’ 

which will result in four coupled equations. In our approach we only consider 
vacuum and thermal fluctuations for the constituent quarks and not for the 
mesons. The explicit expression for the held equations can be found in m 
(Eq. (2)-(5)). 


5. Results and Conclusion 

By solving the held equations Eq. (Q we can investigate the behavior of 
the (pNi (ps order parameters as a function of T at /rs = 0 and as a function 
/xs at T = 0. It was calculated on the lattice [12] that at /xb = 0 the 
value of the pseudocritical temperature Tc should be 151 MeV, while it is a 
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common belief that the order of the transition in at ^ = 0 should be of 
first order. On Fig. [^the and values for all the 40 assignments are 
showrj^ as a function of (the low lying isoscalar mass) together with 
the lattice value for the Tc- It can be seen that in order to be consistent 



Fig. 1. Tc (left) and plb,c (right) values versus ra^L. On the left Tc is shown for 
two different Polyakov potential. 

with the lattice Tc the rurL mass should be below 1 GeV, which can be 
correspond either to /o(500) or to /o(980). However if we would like to have 
first order phase transition on the p-B axis, the mrL mass should be even 

J 0 

smaller (^ 400 MeV). Consequently we investigated the phase boundary for 
parameterizations with relatively small masses. 

On the left panel of Fig. the phase boundary together with the position 
of the critical endpoint (CEP) is shown, while on the right panel the CEP 
variation with the mass is presented. If we increase the mass 
above ps 350 MeV the CEP ceases to exist. 

In conclusion we can say that in order to have a good thermodynamic 
description within the framework of the current model we must have the 
/o(500) particle in the spectrum. 
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Fig. 2. Phase boundary (left) and CEP variation with rrifL mass (right). The CEP 
position on the left is at = 828 MeV, Tcep = 76 MeV) 
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